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Abstract: Electrochemical systems with high capacities demand devices for 
electrochemical impedance spectroscopy (EIS) with ultra-low frequencies (in order of 
mHz), that are almost impossible to accomplish with analogue techniques, but this 
becomes possible by using a computer technique and accompanying digital equipment. 
Recently, an original software and hardware for electrochemical measurements, intended 
for electrochemical systems exhibiting high capacities, such as supercapacitors, has been 
developed. One of the included methods is EIS. In this paper, the method of calculation of 
circuit parameters from an EIS curve is described. The results of testing on a physical 
model of an electrochemical system, constructed of known elements (including a 1 .6 F 
capacitor) in a defined arrangement, proved the validity of the system and the method. 

Keywords: electrochemical measurements; measurement system; electrochemical 
impedance spectroscopy; supercapacitors 



1. Introduction 

Every system may be regarded in a frequency domain displaying frequency logarithm on the X-axis 
and logarithm of module and/or phase angle of transfer function on the Y-axis (Bode plot) [1-3]. For 
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the electrochemical system, the transfer function is, in fact, the impedance of an equivalent electrical 
circuit i.e., its complex form. That is why this method is named electrochemical impedance method 
(EIS) and it is widely used for characterization of electrochemical systems [4—14]. Applying alternate 
voltage, U(S), of different frequencies and constant amplitude to an electrical circuit, responding 
current, I(S), will appear [2,3]. Amplitude and phase angle of this current will depend on voltage and 
impedance, Z(S), of the circuit at that frequency: 

I(S) = 

Z(S) 

where S = ct + jco is Laplace complex variable. For sinusoidal excitation the real part ct is equal to zero, 
so the S becomes S = jco, where frequency, co, is given in s"\ 

Conventional impedance spectra are actually snapshots of transfer functions taken at certain 
well-defined states of the system (usually stationary, constant potential states). However, for a fuller 
description of electrochemical systems the evolution of impedance spectra should be investigated 
during the evolution of the system in both potential and time. With the progress of digital techniques 
this is becoming increasingly possible [15-17]. 

The system used for electrochemical measurements consisted of hardware (PC, AD-DA converter 
NI-621 produced by National Instruments and an analog interface developed at the Technical Faculty 
in Bor) and software for excitation and measurement (LAB VIEW platform and originally developed 
application software) [18,19]. With the goal of achieving a frill mathematical analysis of the measured 
data directly in the Lab VIEW application, it was necessary to develop our own mathematical model 
which is implemented in the measurement software and described in this work. This was not possible 
with some commercially available fitting software packages (EqCwin, Z-view) [20]. 

The possibilities of the software described here are compared with the Thales software of the 
Zahner EIS firm [21,22]. Our model and software are more adequate to the real system because the 
model better describes real electrochemical systems taking into account the complexity of the 
processes. The mathematical model developed herein is adapted to the investigated class of 
electrochemical systems and it is strongly connected with the physical parameters of the system. That 
approach enabled us to obtain analytical values of mutual relationships between the physical 
parameters from the system response and, in that way, to make system optimization following some 
given criteria. This is a significant advantage in compare to the commercial software, where the model 
is not "visible". 



2. Theoretical Part 



By recording amplitude and phase angle of the response current for every frequency value 
(excitation voltage known), one can obtain the module and the phase angle of impedance for that 
frequency; this is presented as one point on the Bode plot which gives the dependence of impedance 
module, Z, on frequency, f, in logarithmic scale. Logarithm is used in a goal to obtain linear 
dependences instead of exponential ones. At frequencies obtained by exfrapolation of straight 
segments, some deviation from straight line appears, and the line slope changes gradually. From the 
heights of the horizontal regions and comer frequencies, one can calculate all the parameters of the 
circuit of which the Bode plot is recorded, i.e. to estimate the equivalent circuit parameters [23-25]. 



Sensors 2009, 9 



7367 



After years of investigating of electrochemical behavior of different electrode materials, different 
equivalent electrical circuits that exhibit the same response on excitations as considered 
electrochemical systems have been found [23-26]. One of the most common was the circuit presented 
in Figure 1 . 

Figure 1. Considered equivalent electrical circuit. 
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Ro corresponds to the resistance of electrolj^e and electrode material, and its value is on an order of 
magnitude of milliohms (mQ) or Ohms (Q). Capacity Co corresponds to double layer formed on the 
electrolyte side. Resistances Ri and R2 (order of magnitude ohm to tens Ohms) are related to slow 
processes of adsorption and diffusion, as well as the capacitances Ci and C2. As a matter of fact, the 
branch RiCi exhibits and describes the inconstancy of parameters in R2C2 branch. R3 is resistance of 
self-discharging, meaning that it is reciprocal to leakage current. Its value is on the order of hundreds 
of Ohms to tens of kiloohms. 

For the adopted equivalent circuit (Figure 1) in a general case the impedance equation is complex 
and not clear enough. So, here a step by step method is applied, one frequency domain after other, 
knowing the nature of the process, i.e. orders of magnitude of the circuit parameters. For very low 
frequencies (on the order of |u,Hz) all capacitors do not conduct electricity, so the impedance of the 
circuit remains the serial connection of Ro and R3: 

Zi = Ro + R3 

where Zi is correlated to the first (the highest) horizontal plateau in Figure 5. At frequencies on the 
order of mHz capacitor C2 conducts, while Ci and Co still are infinite resistances; so, the equivalent 
circuit has the shape presented in Figure 2. 



Figure 2. Equivalent circuit for the second frequency domain (on the order of mHz). 
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The impedance of the circuit presented in Figure 2 is: 

_ S[(Ri +R3) RoC2 +R2R3C2]+Ro +R3 



sCjIRj +R3)+i 



From the conditions for the impedance zero and pole, the comer frequencies may be obtained as: 



27i(R2+R3)C2 



f 2 = 



Rq +R3 



2;e.[(R2+R3)-RoC2+R2R3C2] 
At some higher frequencies (in order of dozens mHz) C2 becomes short circuit, while Co and Ci are 
still in break, so the height of this horizontal region is: 



Z2 = Ro + R23 where R23 = 



R,R 



2 -^^3 



R2 +R3 



At frequencies on the order of hundreds of mHz, Ci starts conducting, Co is still in break, and C2 is 
a short circuit; the equivalent circuit then has the shape given in Figure 3. 



Figure 3. Equivalent circuit for fourth frequency domain (on the order of hundreds of mHz). 
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The impedance of the circuit is then: 

2 ^ S[(Ri +R23)-RoCi +R23R1CJ + R0 +R23 

SC,(R,+R,3)+1 
From the previous equation, comer frequencies may be obtained as: 

f = 1 

2,l.(R,+R^).C, 

and 

' 27i-[(R,+R23)-RoQ+R23R,Cj 

Next horizontal region is obtained at frequencies higher then 1 Hz, when capacitor Ci becomes a 
short circuit, as well as C2, while Co still does not conduct; so it can be written: 
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Z3 = Ro + R123 
where R123 is a parallel connection of Ri, R2 and R3. 

At relatively high frequencies (on the order of kHz) Co starts leading, while Ci and C2 are short 
circuits, so the equivalent circuit becomes as in Figure 4. 



Figure 4. Equivalent circuit for sixth frequency domain (on the order of kHz). 
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The impedance of such circuit is: 



and corner frequencies are: 



2 ^ SRqRi23Cq + Rq + R123 

SR123C0 + 1 



2;iRi23Co 



271R0R123C0 



At the end, the lowest horizontal part of Bode plot is obtained at highest frequencies (on the order 
of tens of kHz) when Co is in short circuit, too, so: 

Z4 = Ro 

The theoretical Bode plot for the whole equivalent circuit given in Figure 1 is presented in Figure 5. 



Figure 5, Theoretical Bode plot for adopted equivalent circuit. 
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3. Experimental 

Testing of the system and developed method was done on a physical model of the electrochemical 
system, constructed of known elements in a defined arrangement as in Figure 1. 

The elements that the physical model was made of were: Ro = 3 Q; Ri = 39 Q, R2 =90 Q; Co = 0,12 
|xF; Ci = 30 mF; C2 = 1,6 F and R3 = 1 kQ. (alternatively R3 = 150 Q). Experiments were performed 
using the following parameters: DC level 10 mV, AC amplitude 5 mV, frequency range 30 \iRz up to 
1 Hz. The obtained curves are presented in Figures 6 and 7. 

Figure 6. Experimentaly obtained Bode plot for the physical model (R3 = 1 kD). 
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Figure 7. Experimentaly obtained Bode plot for the physical model (R3 =150 Q). 
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From the experimentally obtained Bode curve, all parameters of the system have been determined 
by following the next steps: 

From the plateau 4, Ro is obtained immediately from Rq = Z4; 
Horizontal region 1 is equal to Zi, and then R3 can be calculated from: 

R3 = Zi-Ro 

Plateau 2 gives Z2, and then applying: 



R23 = Z2 - Ro and R 2 = ^^^^^ 

R3 +R23 



From horizontal part 3, we get Z3 and calculate R123 = Z3 - Rq. Then Ri can be estimated from: 



_ ^123^23 



' T? -1-1? 

IV23 "I" 1^123 



From the comer frequency fi, capacitance C2 is calculated from: 



C - ^ 



From the corner frequency fs, Ci can be calculated as: 



C 1 
' 27tf3-(R,+R23) 



Finally, from the comer frequency fs, Co is estimated as: 

1 



27lf5-Ri23 



Using the method described above, values of the circuit parameters have been calculated from the 
plot given in Figure 6. The results are compared with those obtained using the commercial software 
EqCwin applied to the data from Figure 6 (Table 1). 



Table 1. Parameters of the investigated equivalent circuit. 



Parameter 


Actual value 


Measured value 


EqCwin value 


Ri [Q] 


39 


41 


45 


Ci [F] 


0.03 


0.03 


0.028 


R2[Q] 


90 


93 


93.4 


C2 [F] 


1.6 


1.58 


1.59 


R3[n] 


1,000 


992 


1,003 
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The plot in Figure 7 gives similar results, except R3, that is, in this case, 150 Q. Plots in Figures 6 
and 7 do not have the fourth plateau for highest frequencies, so Ro could not be determined from such 
a curve. 

4. Conclusions 

Table 1 shows a very good agreement between the actual values of the electrical components 
forming the investigated physical model, the values obtained by the method described in this work and 
the values obtained using a commercial software product. In that way the method, hardware and 
software are fiiUy confirmed. It should be emphasized that this method describes the system very well 
and clearly, but its big disadvantage is that it works with very low frequencies (on the order of |u,Hz), 
that means a need for special equipment (like this described in the present work, or similar) and the 
experiments have a very long duration. The second problem may be resolved by starting the 
experiment from a frequency {2 (much higher than previously indicated), but in that case R3 must be 
determined by some other method (for example potentiostatic). 
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